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In a recent letter [1], making use of quantum Monte Carlo 
(QMC) simulations Rousseau and Denteneer (RD) have in- 
vestigated a two-species Bose-Hubbard (BH) model in a one- 
dimensional optical lattice, with coherent conversion between 
two particles of a species and one of the other, mimicking the 
atom-molecule coherence occurring in cold gases in optical 
lattices close to a Feshbach resonance [21 El. In the notation 
of Ref. 1 the Hamiltonian under investigation is 
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where t a and t m are the atomic and molecular hoppings, U aa 
and U mm are the interatomic and intermolecular interactions, 
U am is the atom-molecule interaction, D is the "detuning" 
between atomic and molecular states (the true detuning being 
D — U aa ) and g is the amplitude of atom-molecule conversion. 
This model was previously investigated by the same authors 
in Ref. 4 for a large variety of parameters. 

One of the main claims of Ref. 1 is the existence of a novel 
phase, dubbed super Mott (SM), which is characterized by 
zero compressibility (namely absence of fluctuations in the to- 
tal particle number N = N a +2N m ) and finite superfiuid frac- 
tion of both the atoms (pj) and the molecules (p™). Moreover, 
the superfiows of atoms and molecules appear to be anticor- 
related, namely the correlated superfiuid density (p c ° r ) turns 
out to be zero. In this comment we argue that the claimed SM 
phase does not contain any superfiuid component, and that the 
individual superfiuid densities of atoms and molecules p a s and 
p™ are not well defined in the model of Eq. ([TJ. In fact, the 
only meaningful superfiuid density is the correlated one p™ r , 
which turns out to be vanishing, revealing a normal phase. 
We corroborate the above statement with the explicit numeri- 
cal calculation of the correlation function associated with the 
coherent counterflow of atoms and molecules, and we show 
that this correlator is short-ranged in the supposed SM region. 

At T = the superfiuid density of one bosonic species is 
defined via the energy cost of a boost in the phase of the op- 
erators associated to that species. E.g. for the atoms, con- 
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^d 2 E /d(f> 2 \j ) =o (where E is the ground state energy and L 
is the size of the system). For a single species BH model, the 
explicit calculation of the superfiuid density within the path- 
integral formalism shows that this quantity can be estimated 
via QMC simulations on a system with periodic boundary 




FIG. 1: Atom-molecule correlator for a L = 100 chain with parame- 
ters: (D/t a , (i/ta) = (3, 8.6) [squares], (7, 9.6) [circles], (13, 10.6) 
[diamonds]. The other parameters are U a — 8t a , U a m — 12i a , 
U m = oo, g — 0.5t a , t m — 0.5t a . According to Fig. 10 of Ref. 4 
all three parameter sets should give a SM phase . 



conditions by means of the fluctuations of the winding number 
(W) of the worldlines associated with the motion of particles 
in imaginary time [5, 6|. When particle-number conservation 
holds — as in the single-species BH model — supercurrents 
are topological invariants (independent of the connected sur- 
face through which the current is measured), and indeed so 
is the winding number for each configuration sampled by the 
Monte Carlo. Yet in the case of Eq. ([TJ, the atom and molecule 
numbers are not conserved separately, so that neither the cur- 
rents of each species are not topological invariants nor the 
atomic and molecular winding numbers, W a and W m . The 
only conserved quantity is the total number N, and hence the 
only meaningful superfiuid density is the one associated with 
the correlated superflow of atoms and molecules, captured by 
the correlated winding number W col - — W a +2W m introduced 
in Refs. 1 4 From a technical point of view, this corresponds 
to the fact that an arbitrary phase boost of the operators, a,j — > 
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, produces a rotation of the atom- 
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molecule conversion term m*-a 2 — > m'-a 2 
grows with the position j (hence leading to an energy varia- 
tion in the ground state which is not infinitesimal in the limit 
4>a,m -^ 0), unless <f> m — 2<fi a , which corresponds to probing 
the correlated response of atoms and molecules, namely the 
correlated superfiuid density. In fact, following Refs.[5]|6]it is 
straightforward to prove that the correlated superfiuid density, 
defined in Ref. .1 as p c s m = \im T ^o{k B TL/2)(W 2 OI ), can 
be written as p c s OT — j^d 2 E / ' d4> a 2 \tj> a =Q with the correlated 
boost (j) m — 2<j) a . Hence the only meaningful superfiuid den- 
sity is the correlated one; which, as shown in Refs. 1 4 is per- 
fectly vanishing in the SM phase. Consequently the supposed 
SM phase is actually normal, and its name is misleading. 



As argued in Refs. Q] |4j the superfluid nature of the 
SM phase stems from counterfiow of pairs of atoms and 
molecules, leading to a zero net mass current. Counter- 
flow superfluidity (CSF) has been theoretically predicted in 
repulsive binary mixtures of bosons (with particle number 
conservation of both species) and it is associated to con- 
densation (or quasi-condensation in one dimension) of com- 
posite objects made of one particle of a given species and 
one hole of the other species (71 El. In analogy to CSF, 
if the SM phase contained counterflowing superfluid com- 
ponents, one would then expect composite objects made of 
two atoms and a molecular hole — or of two atomic holes 
and one molecule — to quasi-condense. We have explicitly 
checked this by investigating the atom-molecule correlation 
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are claimed by RD to correspond to a SM phase. Fig.fTlshows 
the above quantity calculated numerically by a variational Ma- 
trix Product State Ansatz |9| with bond dimension D = 60, 
for a chain of length L = 100 with open boundary conditions; 
all the parameters are chosen so as to be in the SM phase (ac- 
cording to Fig. 10 of Ref. 4) and close to the resonance con- 
dition D = U a . We find that the C am correlator decays expo- 
nentially for all the investigated parameters, as well as all the 
correlators to lower order. Hence no quasi-condensation phe- 
nomenon is observed in the supposed SM phase, confirming 
its fully normal character. 
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